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ABSTRACT 

We examined the morphology-density relations for galaxy samples selected by luminosity and by 
mass in each of five massive X-ray clusters from z = 0.023 to 0.83 for 674 spectroscopically-confirmed 
members. Rest-frame optical colors and visual morphologies were obtained primarily from Hubble 
Space Telescope images. The visual morphologies ensure consistency with the extensive published 
results on galaxy evolution in dense environments. Morphology-density relations (MDR) are derived 
in each cluster from a complete, luminosity-selected sample of 452 galaxies with a magnitude limit 
My < My + 1. The change in the early- type fraction with redshift matches previous work for massive 
clusters of galaxies. We performed a similar analysis, deriving MDRs for complete, mass-selected 
samples of 441 galaxies with a mass-limit of IQ'^^-^Mq. Our mass limit includes faint objects, the 
equivalent of ~ 1 mag below L* for the red cluster galaxies, and encompasses ~70% of the stellar mass 
in cluster galaxies. The MDRs in the mass-selected sample at densities of S > 50 galaxies Mpc~^ 
are similar to those in the luminosity-selected sample but show larger early-type fractions, with a 
weak indication of a shallower slope. However, the trend with redshift in the fraction of elliptical and 
SO galaxies with masses > IQ^^-^Mq differs significantly between the mass- and luminosity-selected 
samples. The clear trend seen in the early-type fraction from z = to z ~ 0.8 is not found in 
mass-selected samples. The early-type galaxy fraction changes much less, and is consistent with being 
constant at 92% ± 4% at E > 500 galaxies Mpc'^ and 83% ± 3% at 50 < E < 500 galaxies Mpc'^ 
Given the mass-limit in our sample, this suggests that galaxies of mass lower than > IO^^'^Mq play 
a significant role in the evolution of the early- type fraction in luminosity-selected samples, i.e., they 
are larger contributors to the luminosity-selected samples at higher redshifts than at low redshift. 
Subject headings: galaxies: clusters: general — galaxies: elliptical and lenticular, cD, — galaxies: 
evolution — galaxies: fundamental parameters — galaxies: photometry 



1. INTRODUCTION 

At low redshift, the cores of clusters of galax ies are 
dominated by mas sive elliptical and SO ga laxies (Smith' 
[1935; Zwicky 1942; Oemleijll97l ). iDresslerl ([i980b) quan- 
tified this as the morphology-density relation (MDR), 
where the fraction of elliptical and SO galaxies (here- 
after early-type galaxies) increases with increasing local 
galaxy density. Much discussion has centered on the as- 
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trophysical reasons for this relation, with contributions 
from many studies. For example, how galaxy morphol- 
ogy depends on th e underlying mass in sta rs has been 
recently explored. iKauffmann et al] (|2003l ) found that 
the concentration of galaxies (a structural parameter cor- 
related with morphology) depends solely on the under- 
lying stellar mass of the galaxies, not on their environ- 
ment. Kauffmann et al. (2004) subsequently found that 
the star-formation rate depends inversely on both the lo- 
cal galaxy density and the galaxy stellar mass. These 
results raise interesting questions about the relative de- 
pendence of galaxy morphologies on galaxy mass and on 
local galaxy environment. 

At hig her redshifts, the star formation rate increases 
rapidly (|Lillv et all [l99fit iMadau et all [1991 . Clus- 
ters of galaxies grow by the accretion of galaxies and 
groups of galaxies, many of which will have a larger star- 
formation rate than galaxies today. Thus, we expect the 
fraction of star-formin g (late-type) galaxi es to increase. 
[Dre s sier et all (I1997D Treu et al. (120031). ISmith et all 



^ 2005V iPostman et al.l (|2005l ). and lDesai et al.l ()2007t ) all 

found an increase in the number of late-type galaxies at 
higher redshifts in clusters and a decreasing fraction of 
SO galaxies at redshifts from the present day out to red- 
shifts z ~ 0.5 — 1. This points to a substantial buildup of 
cluster galaxies in the last half of the age of the universe, 
i.e., over look-back times from 4 to 8 Gyr. 

Endeavoring to understand the details and the causes 
of these changes is challenging, however, since luminos- 
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ity was used to select the galaxies. Variations in the star 
formation rate can cause large changes in the mass-to- 
light ratios of galaxies (these changes can exceed an or- 
der of magnitude) . Galaxies of similar mass thus can be 
spread across a range of luminosities (or be excluded alto- 
gether if too faint) . Selection directly by the stellar mass 
of galaxies would remove some of these problems. For 
example, mass-selected samples would not include very 
low-mass galaxies undergoing bursts of star formation. 
These galaxies would enter luminosity-selected samples, 
and distort comparisons between galaxy populations at 
different epochs. The stellar mass of galaxies changes 
more modestly, even after a major mergers. 

This has particular value for determining the progen- 
itors of the galaxy population in present-day clusters. 
Galaxies of similar masses at higher redshifts are likely 
to form the majority of the progenitors of current galax- 
ies. Except for the rare major merger, mass growth for 
the typical cluster galaxy is likely to be modest. Merging 
and star-formation will move galaxies to higher masses 
over time, so our higher-redshift mass-selected samples 
will not contain all progenitors of z ~ cluster galax- 
ies. However, these processes will impact luminosity- 
selected samples much more dramatically. Proper ac- 
counting and understanding of the evolution of galaxy 
populations (and of the nature of relations such as the 
morphology-density relation) will be better done using 
mass-selected samples. 

Given this, we set out to use mass-selected sam- 
ples to investigate how the morphological mix in clus- 
ters of galaxies changes with local galaxy density and 
how the nix varies with time. The aim of our inves- 
tigation is to determine whether the trends found i n 
the luminosi t y-sele cted samples of iDressler et al.l ()1997D . 
ISmith et all (|2005[ ). and lPostman et all ^20oW appear m 
mass-selected samples. The ultimate goal is to establish 
the progenitors of the z ~ elliptical and SO population 
by examining the masses of z ~ — 1 galaxies. 

We have constructed a sample of five clusters of galax- 
ies; Coma at z = 0.023, CL 1358.4-^6245 (hereafter 
CL 1358) at z = 0.33, MS 2053.7-0449 (hereafter 
MS 2053) at z = 0.59, MS 1054.4-0321 (hereafter 
MS 1054) and RX J0152. 7-1357 (hereafter RX J0152), 
both at z = 0.83. Each cluster has a large cata- 
log of spectroscopic members, combined with high res- 
olution imaging and precise colors in the optical rest- 
frame. These clusters are massive systems with both 
large X-ray luminosities and substantial velocity dis- 
persions. This provides two advantages: first, these 
rich systems provide large samples for studying the dis- 
tribution of morphologies; second, our cluster selec- 
tion minimizes the influence of the relation between 
X-ray luminosity or velocit y dispersion an d early -type 
galaxy frac t ion fo und in IPostman et al. I (I2OOI and 
iDesai et"al] (|2007D . For each cluster, we have rest- 
frame B imaging, generally from the Hubble Space Tele- 
scope (HST), which provides visual morphological clas- 
sificatio ns consist e nt with th ose used in previo u s wor k 
such as lD7essled (Il980bllal): IPostman fcGelleil (I1984D. 
iTreu et all (I2003D: IDr essler et al.l (Il997l) : ISmith et all 
(|2005f) : IPostman et al.l 72005): Desai et al.l (l2007t). Many 
of the c l assifica tions come directly f rom t he work of 
IDressler! (|1980aD and IPostman etall (|2005f ). The un- 
derlying stellar mass is derived using the mass-to- 
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(seelKelson et al.ll2000d: iBefl fc de JomI|2001I : iBefl et al.l 
l2003HHolden et al.ll2006r ) and total magnitudes from our 
measurements. 

We discuss our data in §2, with emphasis on how we se- 
lect galaxies to ensure that we have consistent morpholo- 
gies and photometry across a variety of clusters. We dis- 
cuss our mass estimates as well, with a longer discussion 
in an Appendix IbI In §3, we turn to our estimates of local 
galaxy density, showing that we reproduce the observed 
evolution in the morphology-density relation reported by 
other authors. We then create mass-limited samples of 
galaxies in each cluster and examine the evolution of the 
morphology-density relation for galaxies that, at low red- 
shift, are the massive early- types that dominate cluster 
cores. We discuss our results in §4 with a summary in 
§5. 

Throughout this paper, we assume flm = 0.3, I^a = 0.7 
and Ho = 70 km s~^ Mpc~^. All stellar mass estimates 
are d one using a "diet" Salpeter initial rn ass function 
(IMF iBell fc de JondlMl iBell et al.ll2003l) . This IMF 
has t he slope of a traditional Salpeter IMF (|Salpeteil 
I19551 X = 1.35) over the range of 125-0.35 Mq, with 
a flat slope below 0.35 M©. This IMF produces the 
same colors and luminosity as a Salpeter IMF, with 70%, 
or 0.15 dex, o f t he M / L of the origi nal Salpeter IMF 
(|Salpetedll955D . (|Befl fc de Jong|l200l selected this IMF 
as it matches the properties of low redshift galaxies. This 
I MF yields masses .1—0.15 dex l a rger t han the IMF 
of 'K ennicut^ (|1983f ). iKroupa et all ()1993f ). or IChabrierl 
(20031) . 

2. DATA AND OBJECT SELECTION 

We studied five clusters of galaxies in this paper, and 
they are listed in Table [TJ Each cluster in our sample is 
massive, with a velocity dispersion > 850 km s""'^, and 
is X-ray luminous. All but Coma and RX J0152 were 
included as part of the Einstein Js/Lediwoa Sensitivity Sur- 
vey cluster sample (iGioia et al.lll990l : iHenrv et al.lll992l : 
iGioia fc Luppinolll994D . We list velocity dispersions and 
bolometric X-ray luminosities in Table [2] 

For each of the clusters, we used morphological types, 
photometry and redshifts from the literature to compute 
the fraction of early-type galaxies as a function of lu- 
minosity, mass and local galaxy density. In columns 
(5) — (7) of Table [1] we list the number of galaxies used 
to estimate the local galaxy densities, the number used 
to compute the MDR for luminosity-selected samples, 
and the number of galaxies included in the mass-selected 
samples. When estimating the local galaxy density, 
we selected galaxies down to the limiting magnitude of 
My + 1.5. where My is the absolute magnitude from the 
ISchechten ()1976t ) luminosity func tion in the V band. W e 
used the same value for My as IPostman et al.l (|2005D . 
My = —21.28, and we used the same evolution with red- 
shift where M^ = M^{z = 0)-0.8z. We chose a different 
luminosity- limit, My -\- 1, for our MDRs This limit was 
simply selected to have about the same size as a mass 
selection. We selected a mass limit of 10^'^'^Mq as that 
is our completeness limit at z — 0.83. We explain this 
further in §2.51 Finally, we selected galaxies within 1.25 
Mpc of the cluster center, a limit imposed by the HST 
imaging, except for the Coma Cluster where we extended 
the data out to 1.5 Mpc for comparing our results with 
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previous work in the literature. 

We would like to emphasize that in this paper, early- 
type galaxies are morphologically-classified visually as 
elliptical and SO galaxies from the imaging data. The 
classification was done using the image with the filter 
closest to rest- frame B. Neither the galaxy color nor the 
spectrum was considered, thereby ensuring consistency 
with the extensive body of previous work. We made 
this distinction because strongly bulge-dominated sys- 
tems, such as elliptical and SO galaxies, have been shown 
to evolve differently from spiral galaxies. Furthermore, 
galaxies selected based on spectral energy distribution or 
color often show a mix of morphologica l types and star 
formation rates fsee lPopesso et alll2007l for a recent ex- 
ample) . 

For all the galaxies in our clusters, we collected or de- 
termined rest-frame optical morphologies, B — V colors in 
the redshifted passbands (hereafter B^ — Vz), total mag- 
nitudes and stellar mass estimates. We discuss below the 
data we used to construct these samples. 

2.1. Morphologies 

The morphological classifications came from the liter- 
ature for all the galaxies in our sampl e. For the Coma 
Cluster, we used the morphologies from ' Dressier! (|1980a[ ) 
combined with imaging data from the SDSS. We selected 
early-type galaxies a s those c lassifie d as D, E, E/SO, 
SO/E, SO and SO/a in lDressleil (|1980al ). Late-type galax- 
ies are any galaxies with classifications in the spiral se- 
quence and irregular galaxies. For all of the clusters in 
our sample, we used the same classification criteria. 

For CL 1358, we us ed the morphological types from 
iFabricant et al.l (|200GD . We selected early- type galax- 
ies as those with a classification less than or equal to 0. 
This puts the Sa/SO class in the early-type galaxy cate- 
gory We used the unpublished data of Tran et al. (2003, 
l2005br ) to form a sample of spectroscopic cluster mem- 
bers with morphologies for MS 2053. These classifica- 



tions were done in the same manner as in lFabricant et al.l 
(I2OOOI) and we used the same divisions between early and 
late-type galaxies. Finally, for MS 1054 and RX J0152, 
we used the p reviously p ublished morphologie s from 
iPostman et"alf (I2005D and iBlakeslee et"all (|2006l ). The 
classifications in Postman et al.l (|2005f ) are on a very uni- 
form syst em, one con sisten t with the earlie r samples of 
IFabricant et al.l (|2000( ) and lDressleJ (|1980bl ). 



2.2. Photometry 

We require a consistent set of magnitudes and colors 
to estimate stellar masses . We used the approach de- 
scribed in Blakcsl ee et al.l (|2006l ) to measure color and 
total magnitudes for the cluster galaxies in our sample. 
This p rocedure is a modification of the approach out- 
lined in |^^£okku^^it^n (EMI)- Each galaxy's surface 
brightness distribution was fit by a Sersic model with the 
exponent rangi ng from n = 1 to n = 4, using the program 
galfit (p^ll96llPeng et al.ll2002! ). The fitted model 
is convolved with a point-spread function, accounting for 
the different responses in our data. The fitted magnitude 
was used as our total magnitude estim ate, as we have 
found that to be a reliable measure (Ho lden et all 120041 
2005). We then apphed the CLEAN algorithm (jHogbomI 
1974i) to the original images in all passbands. We mea- 
sured the flux within the half-light radius determined 



from the surface brightness fit to measure the colors. We 
detail below how the colors were measured for each set 
of cluster galaxies since the uniformity and consistency 
with which we carried out the photometry is an impor- 
tant factor in ensuring reliable mass estimates across the 
redshift range studied here. Overall the agreement with 
other studies was very good. 

2.2.1. Coma 

We creat ed a mos aic of SDSS imaging covering all the 
galaxies in iDresslerl ((l980a.) . We fit the Sersic model to 
the g image. For a second estimate of the galaxy colors, 
we froze all of the parameters of the model, except for 
normalization, and fit the surface brightness profiles in 
the u and r data. This is referred to as a model color, 
and is used in the SDSS database. 

We compared the t o tal ma gnitu des with those from 
both Ij^rgensen et al.l (|1992f l and iBeijersbergen et aD 
(|2002allbl l. We found that our total magnitudes were al- 
m ost identical, a d i fferenc e of 0.006 ± 0.017 from those 
of Ei rgensen et al.l (I1992D. However, we found that 
the IBeijersbergen et al.l (2002a. b) magnitudes were 0.12 
mag faint er than our measurements . The r magni- 
tudes from IBeijersbergen et al.l (|2002al lbl) are SExtractor 
MAG_AUTO magnitudes. For the very brightest galax- 
ies, we found that the offset was larger, ~ 0.2 — 0.4 mag, 
which we did not fin d when we comp a red ou r total mag- 
nitudes to those of iJorgensen et all (|1992l ). When we 
measured MAG_AUTO magnitudes from our mosaic, we 
found a similar offset for the brightest galaxies. This 
probably results from the large number of neighbors for 
the brightest galaxies in Coma, suggesting that these 
measurements are not as reliable for total magnitudes. 

We com pared our colors with th ose of the SDSS and 
those from lEisenhardt et al.l (|200 7^ and found very good 
agreement. We found that our colors, measured within 
a half-light radius, were redder by only 0.013 ± 0.004 
mag than the SDSS model colors. We also compared 
our model colors with those of the SDS S, and found that 
those were redder by 0.006±0.005 mag. lEisenhardt et al.) 
(2007) give B — V colors for a number of Coma members. 
We computed a similar conversion that in Table [3l but 
for observed B — V instead of redshifted colors, Bz — 
Vz- Using the relation B -V ^ l-00(g - r) + 0.15 to 
transform our half-light radius g — r colors, we compared 
our calcul ated B — V colors wit h the measured B — V 
colors in Eisenhardt et al.l (|2007f) . We found our colors 
are bluer by -0.008 ± 0.003 with a scatter of 0.029 mag 
when compared with the colors measured within ri, a 
size that roughly corresponds to half of the light, see 
lEisenhardt et all ([2007) for details. 



2.2.2. CL 1358 



For every galaxy in the redshift catalog of lFisher et al.l 
(jl998l), we derived photometry from F606W (Vgoe) and 
F814W (/gu) Wide Fi eld Planetary Camera 2 (W FPC2) 
imaging presented in Ivan Dokkum et all (|1998[ ). We 
used the zero-points from iKelson et alT l 2000a[ r for the 
WFPC2 photometry. We used our standard procedure, 
using the imaging in the /814 passband to determine the 
half-light radii. The colors were m easured in a similar 
manner in Ivan Dokkum et "all (Il998f). We found an offset 
of 0.018 ± 0.006 mag from the Ivan Dokkum et all (|1998D 
colors, and we found a scatter of cr(V606 — ^814) = 0.07 
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TABLE 1 

Summary of Cluster Data 



Observed Rcdshift No. of No. of No. of 

Cluster z Filters Range'' All'' Luminosity'^ Mass'^ 



Coma 


0.023 


u gr 


0.013 <z< 0.033 


202 


100 


95 


CL 1358 


0.328 


Veoe -^814 


0.315 <z< 0.342° 


134 


95 


99 


MS 2053 


0.587 


Veoe -f8i4 


0.57 <z< o.eof 


85 


60 


56 


MS 1054 


0.831 


Veoe *775 2850 


0.80 < 2 < 0.86S 


112 


90 


82 


RX J0152 


0.834 


^'625 *775 Z850 


0.81 <z< 0.87S 


142 


104 


109 



" The redshift range used to define cluster membership. ^ The number of galaxies with redshifts 
in the range of the cluster with My < My + 1.5.^^ The number of galaxies with redshifts in the 
range of the cluster with luminosities My < My + I.'' The number of galaxies with redshifts in 
the range of the cluster with masses > IO'^^Mq.'' The catalog of [Fisher ct al. (1998) contains 
only galaxies in this redshift range. ^ From iTVan et al.l ||200^ ').E From PBlakeslee et al.. (2009) . 



TABLE 2 
Summary of Anoillary Data 



Cluster z km s ^ 10**^ erg s ^ 



Coma 


0.0231 


1008 ± 33=^ 


9.0 ± 0.2t> 


CL 1358 


0.328 




10.2 ± 0.7'^ 


MS 2053 


0.587 


865 ± 71<= 


6.5 ± 0.4^ 


MS 1054 


0.831 


1156 ± 82f 


16.4 ± 0.8S 


RX J0152 


0.834 


919 ± 168^ 


18.6 ± 1.9' 



" [Strublc & Rood (1999)'^ Vocvodkin & Vikhlinin (200l) 
[Fisher et al.i (il998i)'* iHenrvi ()2004 ')° ,Tran et~all 
(1200^) '' ITran et al.l ll2007l ')g IGioia et al.l 

pool )'' IDemarco et~all pOOSi )' IRomer et all pOOOl ) 



mag between the two sets of data. Our total magnitudes 
are -0.021 ± 0.030 mag brighte r than the total magni- 
tudes from lKelson et al.l (|2000a[ ). 

2.2.3. MS 2053 

The photometry we used is ba sed on the archiva l 
WFPC2 Feoe and hiA imaging. IWuvts et al.l (|2004l ) 
gives velocity dispersions and the results of fitting surface 
brightness profiles to the WFPC2 imaging for a subset 
of the galaxies. We measured the apparent magnitudes 
for all cluster galaxies in the WFPC2 imaging using the 
same procedure as for CL 1358. We compared our total 
magn itudes from s urface brightness profiles with those 
from iWuvts et al.l (I2004D. W e find that the /sm magni- 
tudes oflmitsiraOliool) are fainter by 0.037 ± 0.020 
mag than ours, so we consider the two sets of magnitudes 
in good agreement. 

2.2.4. MS 1054 and RX J0152 

iBlakeslee et al.l ()2006l ) measured the colors within the 
half-light radius, with that radius determ ined from mode l 
fits, using the prog ram GALFIT (Pe ng et all I2002D . 
IBlakeslee et alj ()2006( ) let the Sersic exponent range from 
n = 1 to 4. The total magnitud es came from the 
SExt ractor MAG_AUTO parameter (jBertin fc Arnoutd 
|1996[ ). To determine total magnitudes, a 0.2 mag aper- 
ture jx3iTection_was_^^ magnitude 
(sec Holde n et alll2005t IBlakeslee et al.ll2006D . 

.Blakeslee et al.l ( 20061 ) onl y analyzed the then exist- 
ing catalogs of redshifts from IDemarco et al.l (|2005f ) and 



iTran et al.l (|2007t ). We have acquired addition al redshifts 
for bo th clusters since the publication of Blak eslee et al.l 
(|2006f ). We used the Magellan telescopes with both 
the Inamori-Magellan Areal Camera and Spec trograph 
(jBigelow fc Dressleill200l iDressler et alll2006l ) and the 
Low-Dispersion Survey Spectrograph 3. These new 
m embers were ana l yzed i n the same manner as presented 
in Blakesl ee et al.l (20061 ) . 

For MS 1054, we also used the imaging data 
from the Faint Infrared Extr agalactc Survey (FIRES 
iForster Schreiber et all l2006f) . These data cover 
UBVJHK along with the archival WFPC2 photome- 
try for MS 1054. We replaced the WFPC2 data with 
the Advanced Camera for Surveys (ACS) imaging and 
recreated the catalogs. In Appendix El we discuss how 
we processed the data. For most filters, we found very 
small scatter, < 0.01 mag, between our measured mag- 
nitudes and those provided in the FIRES catalogs. Once 
we constructed a photometric catalog, we fit spectral en- 
ergy distributions to measure redshifted colors and the 
mass in stars for cluster members. We used these mass 
and color estimates as a check for our masses and colors 
using only HST imaging. We discuss these in § §2.31 and 
I below. 



2.3. Redshifted Magnitudes and Colors 

To compute the rest-frame magnitudes in Bz and T4 
from the obs erved magnitude s and colors, we used the 
techni que of IBlakeslee et al.l ()2006[ ) and iHolden et al.l 
(|2006f ). We calculated the magnitudes of templates in 
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TABLE 3 

Photometric Conversions 



Cluster 


Redshifted Filter 


Transformation 




Distance Modulus^ 


Coma 




g - 0.65(g - r) + 0.01*' 




35.01 


Coma 


B, - 


1.00(g - r) + 0.15*= 






CL 1358 




/814 -1- 0.27(V606 -/814) + 


0.63 


41.18 


CL 1358 




0.91(y606 -/814)-0.19 






MS 2053 


V, 


/814 - 0.21(V606 - /814) + 


1.07 


42.68 


MS 2053 


Bz - K 


0.70(y606 - ^814) - 0.30 






MS 1054 


V, 


2850 - 0.16(V606 - 2850) + 


0.74b.c 


43.60 


MS 1054 


- Vz 


1.05(1775 - 2850) + O.ll''''^ 






RX J0152 


V. 


2850 - 0.20(r625 - 2850) "t" 


0.78'' 


43.61 


RX J0152 


B, - V. 


1.05(1775 - 2850) + O.ll'' 







The distance modulus assumes = 0.3, Oa = 0.7, Ho = 70 km s~l Mpc~l. ^ The 
Coma gr data and the ACS photometry use an AB zero point. The Veoe filter used in 
this conversion refers to the ACS F606W and not the WFPC2 filter used in CL 1358 and 
MS 2053. 



the rest-frame filters. We then redshifted the templates, 
and computed the magnitudes in the observed filters. For 
the templates, we used exponentially decay ing star for- 
mation rate models from lBruzual fc CharlotI ([2003, here- 
after BC03). These models had exponential time-scales 
of 0.1 to 5 Gyr, covering a range of ages from 0.5 to 12 
Gyr and three metal abundances, 2.5 1 and 0.4 Z©. We 
list the transformations between the observed filters and 
the redshifted ones in Table [3] and include the distance 
modulus to the cluster in our chosen cosmology. For the 
rest-fr ame filters, we used the B and V curves from lBuseil 
(|1978f ) specifically the B3 curve for the B^- Both sets of 
curves are tabulated by BC03. 

For MS 1054 and RX J0152, B^ lies within the ob- 
served wavelength range of the ACS filters, while the 
center of Vz is redshifted to reddest part the zgso filter at 
z = 0.83. To compute the MDR for luminosity-selected 
samples as was done in previous work, we need to com- 
pute the Vz magnitude. We compared our B^ — Vz col- 
ors for MS 1054 with those estimated from the broader 
baseline of photometry from the FIRES survey. We fit 
spectral energy distributions to the FIRES data (see 
Appendix [B]) . We use those fits to interpolate the ob- 
served colors to predict Bz and Vz for each galaxy. We 
found that the Bz and Vz had an average offset be- 
tween the FIRES predictions and the results from Ta- 
ble [3] of 6{Bz — Vz) = —0.006 mag with a scatter of 
<j{Bz — Vz) = 0.005 mag. When we compare the two 
estimates of the Vz magnitude, we find SVz — —0.033 
mag with a scatter of aVz = 0.027 mag. Reassuringly, 
wc find that our redshift colors and magnitudes from the 
ACS imaging alone match what we expect from a much 
broader wavelength coverage. 

We have computed the transformations between the 
observed F606W and F814W WFPC2 photometry and 
the redshifted Bz and Vz magnitudes for MS 2053 and 
CL 1358. These are different conversions from those pre- 
viously listed in the literature. For CL 1358, we found 
good agreement betwee n our transformation s and those 
previously published by Ivan Dokkum et all ^998). We 
found offsets of ~0.1 mag in the colors of red galaxies in 
MS 2053 between our tran sformations , listed in Table El 
and the published ones of lTran et all (|2003D . However, 
our transformation between the observed colors and Vz 



agree with that of IWuvts et all (|2004D to within < 0.05 
mag. For consistency across our data sets, we used the 
conversions we derived that are listed in Table [3l 

2.4. Mass Estimation 

Estimating stellar masses requires a total magnitude 
and an estimate of the mass-to-light ratio (M/L). We 
used the redshifted colors and magnitudes discussed in 
the previous sections to compute the total magnitude 
in the Bz filter, and we used the linear re lation be- 
tween M/Lb and the Bz — Vz color given in iBell et al.l 
(l200l . M/Lb = 1.737(B - V) - 0.942, to derive M/L. 
We assumed the absolute magnit ude of the Sun t o be 
Mb = 5.45." The relati ons of iBell et al.l (I2003D use 
a "diet" Salpeter IMF (see IBell fc de .Tond 120011 . for an 
explanation). This IMF has the same colors and lumi- 
nosity as a Salpeter IMF, but only 70% of the mass. The 
offsets between the derived stellar masses assuming this 
IMF and other pop ular IMF's are given at the end of 
§1 and in IBell et ali (2003) . Generally, most other IMFs 
predict masses 0.1-0.15 dex lower than our estimates. 

The assumption of a single relation between rest-frame 
color and M/L has two likely sources of error. First, 
there is an intri n sic sc atter in M/L at a fixed color, 
which IBell et al.l (120031) estimate as 0.1-0.2 dex. Sec- 
ond, the iBell et al.l l 2003h M/L estimates are appropri- 
ate for z = stellar populations. At higher redshifts, 
star-forming galaxies may have lower stellar masses, as 
they have formed fewer stars , yet they may have simi- 
lar colors. iKassin et al.l ()2006f ) finds some evidence that 
star-forming, field galaxies at higher redshift may have 
lower stellar masses th an galaxies w i th th e same rest- 
frame colors at z = 0. iKassin et al.l (j2006l ) find a shift 
of ~0.1 dex in the stellar mass of star-forming galaxies 
at a fi xed dynamica l mass . This shi ft is large r, ^0.2 
dex, if IKassin et afl (|2006l ) uses the iBell et al.l (|2003l ) 
stellar mass estira ates. For passively evolving galaxies, 
Ivan der Wei et ahl (2006) find that systematic shifts of 
~0.1 dex are possible. To estimate the sizes of the errors 
in our masses, we computed separate M/L estimates for 
a subset of our data. We also compared the M/L from 

This value of Mg varie s slight ly from 
the value in IBinnev fc Merrifieldl l|1998l ) ; see 
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the iBell et all (|2003f ) relations with dynamical mass es- 
timates in all of the clusters in our sample, and we find 
a robust set of mass estimates that are consistent with 
redshift. The results of these tests are given below. 

2.4.1. Comparison of M / L estimates 

For MS 1054, we fit spectral energy distributions to 
the photometric data from the FIRES catalog. In Ap- 
pendix |b1 we discuss our results from fitting spectral en- 
ergy distributions (SEDs), using exponentially declining 
star formation rate, or r, models, to the data of MS 1054 
and a complementary sample of SPSS Fourth Data Re - 
lease (DR4) galaxies HAdelman- McCarthy et"all I2006D . 
For both sets of data, we found t hat the ave r age M /L b 
agrees with the linear relation of iBell et al.l (|2003( ). see 
Figure fTOl We found a scatter of 33%, 0.14 dex, in the 
M/L around the best-fitting relation for MS 1054. Early- 
type galaxies have a very similar scatter of 34%, or 0.15 
dex, while late- type galaxies have 30%. In our lower 
redshift sample created using the SDSS DR4, we find 
a smaller scatter of 28%, or 0.12 dex. These error esti- 
mates agree well with the errors expected from lBell et al.l 
(l2003h . 

2.4.2. Comparison with Dynamical Masses 

We compared our mass estimates with dynam- 
ical masses (Mdyn) from fundan iental plane mea- 
surements from the litera tur e (|j0rgensen fc Franxl 
1991 lJ0rgensen et al.l Il995allhl: iKelson et al.l l2000blla 



Wuvtset al.ll2004HJorgensen et al.ll2005| ). For each clus- 



ter, we combined the published velocity dispersions (u) 
and sizes (reff) for the member galaxies using the rela- 
tion Mdyn = ^Teffa'^/G p0rgensen et al.lll996.) . We find 
that the relations between the dynamical mass, derived 
from the velocity dispersion and re , and the stellar mass 
have consistent zero-points, the difference was ±0.02 dex, 
or ±5% for all of the clusters in our sample (see Fig[TJ) 
The good agreement between the dynamical and stellar 
masses at all redshifts shows that we have a consistent 
mass scale for all early-type galaxies in our sample, and 
one that matches th at of others (van dcr Wcl et al. 2006; 
iRettura et al1l2006f) . 

We do note that the slope between the two mass esti- 
mates appears to be slightly less than one. It seems as if 
the highest mass galaxies have lower than expected stel- 
la r masses. Th i s is a l so seen in some othe r comp arisons 
dCallazzi et al.i [2001 IRettura et all 120061: IKassin et all 
l2006l ). We find th at we recover t he sa me zeropoint for 
this relation as do iGallazzi et al.l ()2006D . after adjusting 
for the different I MF used in th a t pap er, o ver the same 
redshift range as IRettura et al.l ()2006l ) or iKassin et al.l 
(2006( 1. We conclude that there is no systematic shift 
with redshift in our stellar mass estimates, and that our 
photometric stellar mass scale matches the dynamical 
mass scale very well. 

When we compared the stellar mass estimates from 
our colors and total magnitudes to the dynamical esti- 
mates, we found that the scatter increased with redshift 
(see the vertical error bars in Figure [T]) . This increase 
occurs regardless of whether w e compare with the b est- 
fitting relation, the relation of IGallazzi et al.l (|2006l ) or 
a simple slope of one. For Coma, we found a scatter of 
30%, or 0.13 dex, close to the 0.17 dex scatter between 
color-based stellar masses and dynamical ones quoted in 
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Fig. 1. — Photometric stellar masses vs. dynamical masses. 
Red circles are elliptical galaxies, orange squares are SO and SO/a 
galaxies, and blue spirals are later type galaxies. Each cluster has 
a subset of galaxies with velocity dispersion measurements which 
we use to compute masses using the relation Atdyn = ^r^f ftr'^ /G. 
The dotted lines show the approximate mass limits for the velocity 
dispersion samples assuming the magnitude limits and the typi- 
cal colors of the samples. The z = 0.83 panel contains galaxies 
from two different clusters with two different magnitude limits; we 
show only the fainter limit. When computing statistics, we only 
use galaxies with masses above that line. The solid line has a slope 
of one and th e dash ed line shows the best fitting relation from 
IGallazzi et"all ^OOB) using the de Vaucouleurs relation between 
size, velocity dispersion and mass. We show the scatter around the 
relation ofjGallazzi et al. ( 2006) as a vortical error bar, and the av- 
erage error for the dynamical masses as a horizontal error bar. It 
appears that, for the most massive galaxies, masses from the color 
estimates of M/L may systematically, although slightly, underes- 
timate t he dynamical masses. The average for each sample agrees 
with the IGallazzi et all II2006I ) intercept within ±5%. Overall the 
photometric estimates provide a reliable and consistent estimate of 
the stellar mass in galaxies at all redshifts. 

Ivan der Wei et all (|2006l ). At z = 0.83, we found a scat- 
ter of 58 %, or 0.25 dex, clo s e to the scatter found at 
z ~ 1 by Ivan der Wei et al] (|2006[ ) . One possible rea- 
son for the increased scatter is that the errors on the 
velocity dispersion measurements increase. If we remove 
the contribution from the velocity dispersion error in the 
dynamical masses, the scatter only reduces to 55%, or 
0.24 dex, so the increased scatter we find is not from the 
statistical errors on the dynamical mass measurements. 
IKassin et al.l ()2006D find a scatter of 0.16 dex between 
the stellar mass an d the dynamical mass at all redshifts. 
IKassin et al.l ()2006t ) use relations between color and M/L 
that are different from ours, but they recover the same 
slope and zeropoint as we do in Figure [ Tl Th e differ- 
ence in the scatter found by IKassin et all (j2006f ) and our 
results most likely co nies from the details in how the mea- 
surements are made. IKassin et al.l (|2006D explicitly mod- 
els a rotational component for later-type galaxies, which 
they f ound to reduce th e scatter. Both iGallazzi et all 
(|2006[ ) and IKassin et all (|2006D find that the slope of 
the relation between the dynamical mass and the stel- 
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lar mass is less than one . Using that best-fitting slope 
from lGallazzi et all ()2006f ). we reduce the measured scat- 
ter to 25%, or 0.11 dex, for Coma and 46%, or 0.20 dex, 
for our two z = 0.83 cluster s. Such scatter is well in line 
with the expectations from lBell et al.l ([2Q03)- 

2.5. Completeness Corrections and Catalog Limits 

Each data set we used has a different completeness 
limit for the rcdshifts and morphological identifications. 
For each cluster in our sample, we needed to give a weight 
to each galaxy based on the probability that it would 
have either a redshift or a morphological type, based 
on the galaxy's observed magnitude, hereafter w{m,T). 
The same approach was used by Postman et al. (2005). 

For MS 1054, RX J0152, and CL 1358, the spectro- 
scopic redshift catalogs were less complete than our mor- 
phological catalogs. Thus, we weighted the galaxies with 
redshifts based on the fraction of objects with redshifts 
as a function of morpholog ical type and magnitu de, as 
was done in Appendix B of iPostman "eiTan (|200,qD . The 
weight w(m,T), is inverse of the fraction of galaxies 
with redshifts, N{z, to, T) /N{m, T) where N{m, T) is the 
number of objects with classifications and N(z,m,T) 
is the number of objects with classifications and red- 
shifts. Both fractions are calculated in a magnitude bin 
mil > m > mi (we selected bins of 0.5 mag in size). 

The data for MS 2053 have the inverse problem, 
namely complete redshift catalogs, but incomplete clas- 
sifications. Wc therefore weighted each galaxy by the 
fraction of objects classified by the magnitude and 
type, w{m,T), such that w{m,T) = N(mh > m > 
mi,z,T)/N{mh > m > mi,T) where N{m, z,T) is the 
number of late-type or early-type galaxies in a magni- 
tude bin TO/i > m > mi with redshifts and N(z, m) is the 
total number objects in that bin with redshifts. 

We plot, in Figure [2l the stellar mass estimates as a 
function of absolute V magnitude for the five clusters in 
our sample. To create a sample with a uniform-selection 
across all redshifts, we selected a mass limit where we 
were 50% complete at z = 0.83. That limit is M = 
10^ ^-^ Mp, for our ch o ice in IMF (M - W ^- ^^Mp, for 
thelKennicuttI (|1983l ). iKroupa et al.l (|1993l ). or lChabrieil 
(12003D IMF's). 

2.6. Projected Galaxy Density 

To estimate the pro jected galaxy density we used the 
same methodology as jPostman" etall (l2005h . I.e., a rt = 
7 nearest neighbor distance. We included all galaxies 
within 1.5 Mpc of the cluster center with My < My+l.b. 
In total, there are 674 galaxies in all five clusters at or 
above this magnitude limit. 

The local density was calculated to be 

^ " n n2 w„(TO„,r„)~^ 

^ n— 1 

where J^tv is the angular area to the A^th nearest neigh- 
bor. Da is the angular diameter distance for the redshift 
of the cluster, and Wn{mn, Tn) are the weights for galaxy 
n with a given magnitude m and a classification T. We 
discuss how the weights are calculated in ^2.5\ 

Because the weights are used in estimating both our 
densities and our fractions by galaxy type, we derived 
errors for our density estimates by bootstrapping. We 
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Fig. 2. — Photometric stellar masses as a function of absolute V 
magnitude. Red circles are elliptical galaxies, orange squares are 
SO and SO/a galaxies, and blue spirals are later types. The cluster 
samples, are Coma (bottom left), CL 1358 (bottom right), MS 2053 
(top left) and the combination of MS 1054 and RX J0152 (top 
right). The horizontal line is M = 10^^'^' Mq, the completeness 
limit for mass-selected samples. They are derived by determining 
where the completeness if 50% for the reddest galaxies at z = 
0.83 (see i|2.5l for a complete explanation). The vertical line is the 
magnitude limit, Aly + 1, we used for luminosity-selected samples, 
where M* = -21.28 - 0.8^ l|Postman et al.ll2055l ). 

found that the errors on the density estimates ranged 
from 0.17 dex for Coma to 0.24 dex for MS 1054. There- 
fore, when we computed the MDR, we used bin sizes of 
a t least 0.5 d ex. 

iDressled (Il980b[). iDressler et all (Il997t ). and 
IPostman et al.l (|2005D all used a fainter magnitude 
limit. My + 2, than we did when computing the 
densities. We selected a brighter magnitude limit to 
minimize the impact of sample incompleteness. For 
comparison with the se previous results, w e needed to 
correct our densities. IPostman et al.l (j2005l ) had a simi- 
lar problem. Since their data wer e not always complete 
at that magnitude hmit (M^ + 21 . IPostman et all (j2005D 
computed a correction factor for their densities. This 
correction factor was the ratio of the number of galaxies 
in a Schechter luminosity function with a ~ —1.22 down 
to My + 2 over the number of galaxies in the same 
luminosity function down to the magnitude limit of the 
data. Using a similar procedure, and assuming the same 
luminosity function, the densities for a magnitude limit 
of My + 2 would be larger than our de nsities by 0.18 dex. 
In order to co r npare our re sults with iDressled (IT980b) . 
IDressler et all (|1997D . and IPostman et all (|2005f ). we 
have rescaled our densities by 0.18 dex. 

3. THE MORPHOLOGY-DENSITY RELATION AND THE 
EVOLUTION IN THE EARLY- TYPE GALAXY 
FRACTION 

We combined the data discussed above and created two 
samples for each cluster in order to derive MDRs and 
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to study the evolution in the early-type galaxy fraction 
with redshift. First, we replicated the previous work us- 
ing luminosity-selected samples, then we selected galax- 
ies using the masses estimated above. For each sample, 
we estimated the MDR in a consistent manner, using 
the same sets of densities for both luminosity- and mass- 
selected samples. In the following subsections, we outline 
how we derived MDRs for our clusters. 

3.1. Early-type Fraction Estimates 

The early-type fractions for the clusters were not sim- 
ply based on an integer number of objects. Instead, the 
early-type fraction was the sum of weights of elliptical 
and SO galaxies above our magnitude limit or mass limit, 
divided by the sum of the weights of all galaxies meet- 
ing that same limit. In general, these weights are close 
to one, and not using these weights changes the result- 
ing fractions by only a few percent. While the changes 
are small we felt it important to be consistent and thor- 
ough in our approach. The errors on the fractions range 
from 3% to 15%, depending on the number of galaxies. 
These errors are close to, but slightly larger, than those 
expected from using Poisson errors. We added in quadra- 
ture an error of 6% on the early-type fractions we find 
for RX J0152 a nd MS 1054. This is t he rms classification 
error found bv IPostman et all ()2OO50 . 

3.2. Luminosity-Selected Samples 

Our first step was to create luminosity-selected sam- 
ples to compare with previous work on the MDR. As 
discussed in fJJJ we computed the morphology-density 
relation down to a magnitude limit of M *^ -\- 1 for all 
cluste rs, the same limit as was selected by ISmith et al.l 
(|2005f ). We plot the resulting morphology-density rela- 
tions for a sample of a total of 452 galaxies in Figure 
[3l We tabulate the plotted fractions in Tabled We see 
MDRs in our clusters. Our data clearly shows a decrease 
in the early-type fraction with decreasing density for 
luminosity-selected samples. This evolution is consistent 
with t h at shown in previou s wo rk such as Dressie r et al.l 
(1991 . ISmith et all (|2005h and lPostman~et all (I2005h . 

3.2.1. Evolution of the Luminosity-Selected Early-type 
Fraction 

The MDRs in Figure |3| also show a clear trend of 
decreasing early- type fraction with redshift. To quan- 
tify this we examined the fraction of early-type galaxies 
in two density bins as a function of redshift, S > 500 
galaxies Mpc~^ representing high galaxy densities, and 
500galaxiesMpc~^ > S > 50 galaxies Mpc~^ for low den- 
sities. The cutoff of S = 500 galaxies Mpc~^ was chosen 
to split the combined MS 1054 and RX J0152 into two 
equal-sized subsamples. The lower limit is the density 
that contains 95% of our sample. We plot our early-type 
fractions in each of these density bins as a function of 
cluster redshift in Figure [4l We find a trend of a de- 
creasing early-type fraction with redshif t . as seen before 
(IDressler et al.ril997t ISmith et all 120051 : jPostman et all 
|2005() . For comparison, we show the average fractions 
in the same den sity bins from the low-redshi f t sam ple of 
Dressled (ll980bD. up dated in IDressler et all ()1997[ ). the 
Dressier et al.l (Il997f) z ~ 0.4 sample and the z ~ 1 sam- 
ple of IPostman et al.l ()2005l ). Our data agree with the 
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Fig. 3. — Morphology-density relation for luminosity-selected 
galaxies by cluster. We selected all galaxies to an absolute V mag- 
nitude of M^-l-1 or Mv = -20.28-0.82, (see Postman et al. 2005), 
for a total of 452 galaxies in all five clusters. The cirlces connected 
by lines represent Coma (blue circles), CL 1358 (green circles), 
MS 2053 (orange circles) and the combined sample of MS 1054 
and RX J0152 (red circles). The bins are shown along the bot- 
tom, with triangles marking the centers. The galaxy fractions are 
derived in the same density bins, but are spread out for clarity 
around the bin centers in redshift order. The data for Coma con- 
sist of the morphologies from Dressier ( 1980a) but photometry from 
the SDSS. For comparison, the updated Dressier (1980b) results - 
from Dressier et al. ( 19971) - are shown as squares, the results of 
jPresslcr et al. (1997) at z ~ 0.4 are shown as triangles, and the 
IPo stman et al. (2005) results at z ~ 1 are shown asblack circles 
with error bars. We find good agreement between our morphology- 
density relations and the previous results from the literature. A 
trend in the early- type fraction with redshift is apparent, and is 
quantified further in Figure |4] We list the values for this figure in 
Tabled 

trend shown in the other published work, given the size 
of our sample. 

The dashed line in Figure [D is a fit to the data that 
includes errors on the estimated early-type fraction. We 
find slopes that differ from zero by more than 2a for both 
fits; i.e., the non-zero slope is significant at the 95% confi- 
dence limit. At high densities, we find that the early-type 
galaxy fraction decreases by -0.20z ± 0.10 z, evolution 
at the 95% confidence limit. For the low-density sample 
we find -0.292: ± 0.10 z, evidence at the 99% confidence 
limit for evolution. Generally, studies find that most 
of the evolution at low to moderate redshifts ha ppens 
at lower densitie s (|Treu et al.l l2003t ISmith et al.l 120051 : 
iDesai et"aI]|2007D . a result we replicate with our data. 
Our results at all densities are clearly consistent with th e 
trend seen from the prev ious work of l^issle^ (|1980bf) . 
IDressler et all (|1997f ) and IPostman et al.l (|2005D . We list 
the fractions in the columns (3) and (4) of Table jS) 

3.3. Mass-Selected Samples 
3.3.1. The Mass-Selected MDR 
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Fig. 4. — Early- type fraction as a function of redshift in two 
density bins for luminosity-selected samples. The density bins are 
listed on the figure in galaxies Mpc~^ (and differ from the bins used 
in Figure [Sjl. Blue, green, orange, and red circles represent Coma, 
CL 1358, MS 2053 and the two z = 0.83 clusters, respectively. For 
comparison, the squares at lo w redshift are the ave rage fractions 
in the listed density bin from IDressler et al.l I I1997I) . an updated 
version of [Dr essler ( 1 980bl), the t riangl es show the average fraction 
at z ~ 0.4 from IDressler et al.l 1 119971 ) , and the black circle s rep- 
resent the high-redshift fractions from [Postman et al.l II2005I ). We 
note that these are not weighted by the number of galaxies at each 
density, but are simply the average early-type fraction inside the 
range of densities listed in the panel. The dashed line represents 
a fit to the data from this paper, see Table \E\ We see a strong 
trend in the evolution of the early- type fraction with redshift, and 
a lower early-type fraction at lower densities at most redshifts. The 
evolution with redshift in the early-type fraction is significant at 
the > 95% confidence limit and is similar to that see by others. 

The next step was to determine the MDR for mass- 
selected samples. We used the same weights and densities 
as in Figure [3] and H3.2I Our redshifts or morphological 
classifications became incomplete at M < 10^°'^ Mq, so 
we select this as our mass limit. There are a total of 441 
galaxies above that mass limit in all five clusters. This 
mass limit corresponds to a factor of 4 belo w z = M.* 
for high-density regions (Ba ldrv et al.ll2006h ^^. assuming 
no evolution in the mass function. If we integrate over 
this mass function, our mass limit encompasses 68%-73% 
of the total stellar mass of the cluster population, de- 
pending on the mass limit for the galaxy population (10^ 
or IO^'^Mq, respectively). To c ompute this, we used the 
best-fitting mass function from iBaldrv et al.l (|200(l ) for 
high-density regions, i.e., S > 10 galaxies Mpc~^. Thus 

As noted above, we use an IMF that adds 0.15 dex to the 
masses of galaxies, raising the value of log A^* from 11.05 to 11.2. 



CO 
Sh 

CD 
ft 



CO 




50 100 200 



500 



S (gal Mpc ^) 



1000 



Fig. 5. — Morphology-density relations for mass-selected samples 
of galaxies in our clusters. The sample consists of all galaxies in 
our five clusters ab ove a mass limit of > 10^°-^Mb (> 10^° "^^ 
for th e IMF's from IKennicuttI IT983: IKroupa et al.lll993l: IChabrieil 
I2003I V or 0.25 M* for a total of 441 galaxies in all five clusters. Our 
mass-selected sample thus extends significantly below Ai* (and 
L* , since the mass limit corresponds to 1 mag below L* for the 
reddest cluster galaxies), and includes ~ 70% of the stellar mass 
in galaxies in the cluster. Th symbols are the same as in Fig. |3] 
The samples of galaxies at E < 50 galaxies Mpc~^ are too small to 
be statistically significant. In contrast with Fig. |3]the early- type 
fraction is larger, and the morphology-density relation is less clear. 
The trend with redshift also is noticeably less; see Figure [3 We 
list the values for this figure in Table |5] 

our mass-limited sample encompasses a wide range of 
masses for cluster galaxies, and includes the majority of 
the stellar mass in galaxies in clusters. 

We show the resulting mass-selected MDR in Figure [5l 
The mass-selected MDR (Figure [SJ is different from the 
luminosity-selected MDR (Figure [3]). First, the slopes 
appear to be slightly less. Second, at low redshift the 
MDRs appear roughly comparable, but at higher red- 
shift, at a given density, the mass-selected early- type 
fraction is 10%-20% higher than the luminosity-selected 
fraction. The contrast is quite striking, considering that 
our mass sample extends significantly below J^* (and 
below L* , since the mass limit corresponds to ~ 1 mag 
below L* for the reddest cluster galaxies). See Tables [J] 
and[5]for the resulting fractions in comparison with those 
from luminosity-selected samples. 

The MDR we show in Figure[5]is somewhat of a hybrid. 
We estimated the galaxy densities using magnitude- 
limited samples, while the early-type fractions came 
from mass-limited samples. We used a luminosity- 
selected sample to estimate the local galaxy to minimize 
the changes in comparing the two samples (luminosity- 
selected and mass-selected) . In Figure [51 we show the 
morphological fractions as a function of a galaxy density 
estimated using a mass-selected sample. The mass limit 
we used was IO^'^'^Mq and we used the same weights 
as discussed in i|2.5l We find that the mass-selected lo- 
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Fig. 6. — Morphology-density relation for mass-selected samples 
of galaxies using mass-selected galaxy densities. The symbols are 
hte same as in Fig. \3\ This figure is similar to Figure |5] where 
we plot the fraction of elliptical and SO galaxies above the mass 
limit of > 1O^O-^M0. In this fi gure, the local galaxy densities are 
computed from the mass-selected sample (Em), not the luminosity- 
selected sample as in Figure \5\ The change in our estimate of the 
local galaxy density does not change our results. The typical values 



luminosity-limited sample, which lvan der Wei et all 1)20071 ) finds as 
well for field galaxies. 

cal galaxy densities are generally a factor of ~2 lower 
than those from luminosity-sele cted samples for the same 
galaxy, similar to that found bv lvan der Wei et aP (l2007h 
for field galaxy samples. This offset in the density does 
not, as expected, change the overall shallow slope, we 
find the same value for the slope as we do in Figure [3 
When plotting against these new densities, we still find 
a lack of evolution as we found in Figure [5l 

The mass-selected MDR appears more shallow than 
the luminosity-selected relation. We tested this in a sim- 
ple way. We fit the data from Dressier ( 198Gb) with a 
linear slope in the density range our sample covers. We 
compared that slope to the average slope for the MDR 
of all five clusters in our sample that we plot in Figures [5] 
and [5] The measured slopes are 0.26 ±0.09 forlDressleif 
(jl98Gb) and 0.25 ± 0.13 for Pos tman et all (|2005l ). Our 
mean luminosity-selected MDR has a slope of 0.18±0.07, 
while our mean mass-selected MDR have a mean slope of 
0.12±0.08. To measure these slopes, we average the sam- 
ples after removing the offset between the average early- 
type fractions for each cluster. This preserves the shape 
of the MDR for each sample. The slope in the mass- 
selected MDR appears to be less than our luminosity- 
selected MDR, but it is not statistically-significant. We 
find a larger, but still not statistically-robust difference 
bet ween the slope of the r elation from Dressier (1980b) 
and iPostman et all ()2005l ). and our mass-selected sam- 
ple. 

Overall, the mass-selected MDRs do not differ substan- 
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Fig. 7. — Early-type fraction as a function of redshift for mass- 
selected samples in two density bins. The symbols are the same as 
Fig. |4] The density range in galaxies Mpc~^ and the mass limit 
are labeled in each plot. We show the average early-type fraction 
in each bin as a solid line, with the height of the bar showing 
the error on the average. We find mean early-type fractions of 
83 ± 3% at 500 > S > 50 galaxies Mpc-^ and 92 ± 4% at 
S > 500 galaxies Mpc~^. The dashed line is the best-fitting 
linear relation to the change in fraction with redshift, as in Fig. 
U] We list the fractions in Table [6] There is little evidence for 
evolution in the morphology-density relation above the mass limit, 
with the best-fitting slopes half or less than half those of Figure 
U] The evolution in the luminosity-selected samples must occur 
predominately in galaxies with M < 10^'^'^Mq. 

tially from the luminosity-selected MDRs, except in one 
very important aspect - their evolution with redshift. 

3.3.2. Evolution of the Mass-Selected Early-type Fraction 

Again, as for the luminosity-selected sample, we split 
the sample into high- and low-density bins to establish 
the trend with redshift in the early-type galaxy frac- 
tion. We measured the mass-selected early-type frac- 
tions in the same two density bins as Figure [H This 
is shown in Figure [T] Clearly the trend with redshift 
seen in Figure [4] is much weaker in Figure [T) At high 
densities, we find that the early-type galaxy fraction de- 
creases by -0.13 ± 0.09 z while for the low density sam- 
ple we find -0.14 ± 0.10 z. These slopes are shown as 
dashed lines in Figure [71 Neither represents evolution 
at the 95% confidence limit. We note here that the 
errors are the same size as those given earlier for the 
luminosity-selected samples, but the slopes are half the 
size. We find mean early-type fractions of 83 ± 3% at 
500 galaxies Mpc~^ > S > 50 galaxies Mpc"^ and 92 
± 4% at S > 500 galaxies Mpc~^. This reaffirms the 
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TABLE 4 

Luminosity-Selected Morphology-Density Relations 



Cluster 



1580 > S > 500=^ 
galaxies Mpc~^ 
(%) 



500 > S > 167=' 
galaxies Mpc~-^ 
{%) 



167 > S > 50=^ 
galaxies Mpc~^ 
(%) 



Coma 
CL 1358 
MS 2053 
MS 1054 



RX J0152 



0.023 95 ± 5 

0.328 84 ± 9 

0.587 74 ± 12 

0.83 83 ± 8 



95 ± 5 
82 ± 6 
65 ± 10 
69 ± 8 



75 dz 8 
74 ± 9 
47 ± 12 
69 ± 11 



The fraction of elliptical and SO galaxies with luminosities My 
range. 



<M*- 



■ 1 in the density 



TABLE 5 

Mass-Selected Morphology'-Density' Relations 



Cluster 



1580 > S > 500=^ 
galaxies Mpc~^ 
(%) 



500 > S > 167=^ 
galaxies Mpc~^ 
{%) 



167 > S > 50=^ 
galaxies Mpc~^ 
(%) 



Coma 0.023 96 ± 5 99 ± 3 80 ± 8 

CL 1358 0.328 87 ± 7 88 ± 5 79 ± 8 

MS 2053 0.587 77 ± 10 81 ± 9 60 ± 13 

MS 1054 & RX J0152 0.83 88 ± 7 75 ± 8 82 ± 10 



The fraction of elliptical and SO galaxies with masses > IO^'^'^Mq in the density range. 



lack of evolution seen in Figure [5] and Figure [SI We show 
the mean value of the early-type fraction as a solid line. 
The early-type fractions for the mass-selected sample in 
Figure [7] are ~20% higher than the z > 0.3 results from 
the literature shown in Figure [H and ~10-20% higher 
than our own results from the luminosity-selected sam- 
ples. These results can be seen in fifth and sixth columns 
of Table El as well. 

We measured the slopes for the mass-selected local 
galaxy densities in two density bins as well. We sim- 
ply use 250 > S > 25 galaxies Mpc~^ and E > 
250 galaxies Mpc^^ as the typical mass-selected density 
is half that of the luminosity-selected density. We find 
slopes of -0.13 ± 0.12 z at high densities and -0.16 ± 
0.12 z at low densities. The fraction of early- type galax- 
ies in the high density is 89 ± 4% while we find 82 ± 3% 
in the low-density bin. These are the same within the 
errors as was found for the early-type fractions in the 
luminosity-selected density bins. 

The lack of a trend in the mass-selected sample points 
to galaxies of mass less than 10^^-^ Mq as being the likely 
contributors to the evolution at high redshift seen in the 
luminosity-selected samples, as shown in Figure [3] and 
Figured] This would be consistent with expectations from 
the increasing global star formation rate at higher red- 
shift. The resulting substantial changes in their M/L 
would lead to low-mass galaxies brightening into the 
luminosity-selected samples. 

3.3.3. Is Coma Unusual? 

A minor but interesting issue arises regarding Coma 
and its somewhat higher early-type fraction. In Figures 
[3 and [51 Coma appears to have a larger number of 
early-type galaxies at the highest densities. We find that 
Coma has an early-type fraction 9% ± 6% higher than 
the average early-type fraction in the highest-density bin 



in Figure [7| The apparent higher fraction of early-type 
galaxiess and the lack of any late-type galaxies at very 
high densities for Coma could have several causes. First, 
the number of spiral and late- type gala xies in th e other 
clusters could be overestimated. Diafe rio et al.l ()200lD 
provides a cautionary note, stating that in some sim- 
ulations, star-forming field galaxies could be projected 
into a cluster even for redsh ift-selected samples. Us - 
ing the data from GOODS (jvan der Wei et all [20071 ) . 
we estimate that 2±2 late-type galaxies above our mass 
threshold per cluster would fall within the redshift win- 
dows listed in Table [T] This would change the early- 
type fraction by 2% (for MS 1054 or RX J0152) 10% 
(for MS 2053), depending on the richness of the cluster 
and the density bin in question. This would raise the 
higher redshift cluster early-type fractions, to an aver- 
age of 89%, closer to the Coma value. A second, more 
interesting possibility, is that there is evolution in the 
number of spirals in the cores of clusters from z ~ 0.3 to 
z ~ 0. This would be mild evolution, 9 ± 6% at the high- 
est densities we sample, as compared with the changes of 
25%-30 % seen in prev i ous lu r ninosity-selected samples, 
such as [Dressier et all ([1997[ ). [Postman et al.l ([2OO50 or 
[Desai et al.l ([2007[ ). Third, Coma could be on the high 
side of the distribution of early-type fraction among low- 
redshift clusters, as indicated in Figure [31 The signifi- 
cance of contributions from these effects could be eval- 
uated better with larger numbers of clusters and larger 
mass-selected galaxy samples. 

3.3.4. The Mass Dependence of the Early-type Fraction 

As a further check on the lack of evolution, we in- 
vestigated whether, in some particular mass range, we 
could see evolution in the early-type fraction. We subdi- 
vided our sample into different mass bins, and derived the 
early-type fraction in Figure [8l We found no significant 
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TABLE 6 

Early-type Galaxy' Fractions in Two Density Bins 



Cluster E > 500=^ 500 > S > SO'' S > 500*^ 500 > S > SO"" 

galaxies Mpc~'^ galaxies Mpc"'^ galaxies Mpc~'^ galaxies Mpc"'^ 
(%) (%) (%) (%) 



Coma 


0.023 


96 ± 4% 


86 ± 4% 


96 ± 4% 


86 ± 4% 


CL 1358 


0.328 


86 ± 8% 


80 ± 5% 


89 ± 6% 


85 ± 4% 


MS 2053 


0.587 


74 ± 12% 


57 ± 8% 


77 ± 11% 


74 ± 7% 


MS 1054 & 


i RX J0152 0.83 


83 ± 7% 


69 ± 8% 


89 ± 7% 


77 ± 7 



^ The fraction of elliptical and SO galaxies with luminosities My < My + 1 in the density range.'' The 



fraction of elliptical and SO galaxies with masses > in the density range. 



TABLE 7 

Early- type Galaxy Fractions in Mass Bins 



M > 10"-2 1O"''M0 lO^^-^M© lO^'-'^M© 
Cluster z % % % % 



Coma 0.023 

CL 1358 0.328 

MS 2053 0.587 

MS 1054 0.83 
and RX J0152 



89 ± 6 96 ± 4 

86 ± 10 100 ± 5 

79 ± 14 80 ± 14 

85 ± 9 87 ± 8 



88 ± 7 89 ± 5 

75 ± 8 85 ± 6 
75± 9 71 ± 13 

76 ± 9 87 ± 8 



^ The fraction of elliptical and SO galaxies in the listed mass range at all local galaxy 
densities; the binwidth is 0.2 dex except for the highest mass bin, which includes all 
galaxies above the mass limit. 




Mass (Mq) 



Fig. 8. — Morphology-mass relation for all five clusters. Each 
mass bin is 0.2 dex in width, starting at lO'" '' Mq. The data 
points represent the fraction of early-type galaxies with Coma (blue 
circles), CL 1358 (green circles), MS 2053 (orange circles) and both 
MS 1054 and CL 1358 (red circles). The data points are offset in 
mass about the mean in the bin for clarity. The triangles show 
the bin centers. The last bin contains all galaxies with masses 
> lO'^'^ ^Q. We examined this distribution to see if the early- 
type galaxy fraction did change at, for example, the lower end of 
our mass-distribution. We find no statistically significant slope in 
the early-type galaxy fraction with mass, and the average fractions 
agree with those we find in different density bins, see Table [7] The 
constant early-type fraction we see in Fig. \5\ Fig. \E\ and Fig. [7] 
comes from contributions over the whole of the mass function, and 
is not caused by just the peak of the mass function is dominated 
by early- type galaxies. The typical mass of a late- type galaxy 
contributing to Figures \3\ and |3] must be below our mass limit of 
10^°-^ Mq. 

change of the early- type fraction with mass. Linear fits to 
the data in Table [71 the fraction of early- type galaxies in 
each mass bin, yield slopes of 0.09 ± 0.12 for Coma, 0.23 
± 0.15 for CL 1358, 0.15 ± 0.30 for MS 2053, and 0.01 ± 
0.13 for the combined sample of MS 1054 and RX J0152. 



For all five clusters, we find no statistically-significant 
variation from a slope of zero. The offsets in the mean 
with increasing cluster redshift can also be seen to be 
very small (consistent with the lack of evolution seen in 
Figure [7] for the mass-selected sample). In general, the 
lack of a change in the early-type fraction seen in Figures 
[Hand [7] appears to be true for the whole of the sample, 
not just a subset of galaxies in a certain mass range. 

We also tested whether changing the mass limit at 
lower redshifts made any difference. This arises from a 
potential concern that we may have introduced a system- 
atic bias in our results by keeping a constant mass limit 
for all redshifts, given the expectation of some system- 
atic increase in the mean mass of galaxies (from merg- 
ing and/or star formation). We increased the limit for 
Coma by a factor 2 x to correspond to a significant mass 
buildup between z ~ 0.8 and z ~ 0. This would ac- 
com modate any likely mean mass buildup (e.g. for the 
field iBell et alJlMft finds that - 50%-70% of today's 
galaxies with masses > lO^^'^M© may have undergone 
mergers since z ~ 0.7). The changes in the early- type 
fraction were insignificant, as would be expected from 
Figure [H This is discussed further in §4.1. 

4. RESULTS AND DISCUSSION 

4.1. Mass-Selected and Luminosity- Selected Samples 

As discussed in the §1, we used mass-selected samples 
because we expect only modest and systematic evolution 
in galaxy masses with time, unlike luminosity-selected 
samples where galaxies can enter and leave the sample 
as bursts of star formation dramatically increase the lu- 
minosity for a period. A mass-selection provides a more 
robust means of identifying likely progenitors of z = 
cluster galaxies in higher redshift clusters. However, 
even the mass-selected sample of higher redshift cluster 
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galaxies will not contain all the progenitors for two rea- 
sons. First, galaxies will continue to grow in mass, either 
through star formation or by merging. Mass growth will 
bring objects from below our mass limit at high redshift 
into the sample at low redshift. Second, the ear ly-type 
fraction is lower in the surrou nding field (.Noeske et al.l 
I2007t Ivan der Wei et al.ll2007l ). Field galaxies will con- 
tinue to infall into clusters and build up the cluster pop- 
ulation. Thus, even if there is no evolution from z — 0.83 
to the current epoch in the fraction of early-type galaxies 
in mass-selected samples, we do not expect that the total 
mass in early-type galaxies remains constant. The total 
mass in early-type galaxies in the cluster environment is 
likely to increase (as has bee i i suggested for the field by 
iBell et al.l[200l iFabeil [20071 : iBrown eFall [20071 for ex- 
ample) - it just has to increase through processes which 
maintain a similar high fraction of early-type galaxies at 
masses > lO^^'^M© at any time over the last 7 Gyr. 

It is interesting to explore the types of galaxies in our 
samples. Our mass-selected samples of cluster galaxies 
show little or no change in the early-type galaxy frac- 
tion with redshift for masses > IO'^'^'^Mq, i.e., for masses 
above ~25% of the z = characteristic mass, or Mi,. 
This mass limit encompasses the majority (~70%) of 
the stellar mass in galaxies in these clusters. When we 
use the same parent sample but select galaxies by lu- 
minosity, we recover the MDR and the evolution seen 
by other authors. We illustrate this selection for the 
combined sample of MS 1054 and RX J0152 in Figure 
[9| where we show the distribution of masses and galaxy 
morphological types over a range of densities for galax- 
ies in a luminosity-selected sample (all galaxies have 
My < + 1). At masses lower than 10^°'^ Mq, 
this sample of galaxies becomes dominated by late-type 
galaxies. These lower-mass, late-type galaxies do not fall 
in our primary sample - because we are incomplete at 
these low masses - but do fall into the luminosity-selected 
sample, and so must be the primary contributers to evo- 
lution in luminosity-selected samples. 

4.2. Merging and Mass Assembly 

Merging may play an important role in the forma- 
tion of massive, red gala xies between z ~ 1 and today. 
Ivan Dokkum et al] (|1999D found that the merger rate in 
MS 1054 was significantly higher than in z ~ clusters, 
and that the l a rge m ajority of these are "dry" merg- 
ers. iBell et all ((20061 ) finds that, in the field, as many 
as 50%- 70% of today's galaxies with masses > lO^^-^M© 
could have un d ergon e mergers since z ~ 0.7, although 
iMasiedi et al.l (|2006f) finds little evidence of mergers 
among massive, non-star forming galaxies at z < 0.36. 

Merging raises an interesting issue. Could a signifi- 
cant fraction of the z ~ cluster population lie below 
our mass limit at high redshift? As the following argu- 
ments show, this appears to be unlikely. One estimate of 
the merger rate is the n umber of close pairs o f galaxies. 
iTran et all (|2005aD and [Postman et~aLr (|2007D searched 
for close pairs of galaxi es in MS 1054 and R X J0152. 
Close pairs are defined in lPostman et al.l (j2007[ ) as galax- 
ies within 30 h~^ kpc, but more than 200 h~^ kpc from 
the brightest cluster ga laxy. Our samp l e incl udes a sub- 
set of those pairs from iPostman et al.l (|2007f ). as we re- 
quired both galaxies to be members of RX J0152, and 
the galaxies to be separated by less than 300 km s~^ in 



the rest-frame of the clust er, the same criteri on of the 
whole sample discussed in iTran et all (|2005af ). These 
two samples of potential mergers, 14 in MS 1054 and 9 in 
RX J0152, resulted in 23 galaxies in pairs or triplets likely 
to merge in the next Gyr, or 12% of our My < My + 1. 
sample. The fraction of mer ger candidates in the res t 
of our sample is small, <1% ([van Dokkum et allll999f) . 
pointing to the possibility that cluster galaxy mergers 
were more prevalent either at higher redshift or in cer- 
tain kinds of clusters. Nonetheless, 12% could represent 
the rate of merging in z ~ 0.8 clusters. 

We show the 23 members of these close pairs of galax- 
ies as black open circles in Figure [9| The average mass of 
a galaxy in these pairs, for galaxies above our luminos- 
ity limit, is 10^^-^^^-^Mq. The mean early-type galaxy in 
the whole sample is the same mas s, 10^^-^^°-^Mfn. T here- 
fore, the sampl e of clo se pairs in [Tran et all ([2005al) and 
IPostman et al.l ([2007D appear to have masses typical of 
Mir at z ~ cluster early-type galaxies. It is possible 
that the progenitors of some low-mass (close to our limit) 
z ~ early- type galaxies have a mass at z ~ 0.8 such 
that they would faU below our > lO^'^ '^Mo hmit. The 
close pairs we observe suggest this possibility, but most 
(19 of 23) of the galaxies observed have masses signifi- 
cantly more than a factor of 2 above our mass threshold. 
Second, only two are spirals, pointing to a prevalence of 
"dry" merging. Assuming no new galaxies, "dry" merg- 
ing will decrease the number of early-type galaxies above 
our mass-limit. However, galaxies below our mass limit 
could also merge. The overall impact from merging on 
the early-type galaxy population depends on a number of 
factors, and requires up-to-date galaxy formation models 
to disentangle. However, if the close pairs in our sample 
are representative of the typical mergers between mas- 
sive galaxies, the progenitors of the z ~ galaxies near 
our mass limit may not have built up through massive 
mergers and merging will not have dramatically changed 
the early-type galaxy fraction. 

As mentioned in §3 we carried out a test to see if setting 
the mass limit in Coma at a higher value, > IO^^'^Mq, 
made any difference (compared to > IO^^'^Mq for the 
z — 0.83 clusters). This would correspond to a large o:^ 
2times mass buildup from z ~ 0.8 to ~ (c.f. [Bell et aD 
[20061 results mentioned above). Not unexpectedly, the 
results are unchanged (see Figure [H|). The early- type 
galaxy fraction increases, as would be expected, but very 
modestly by 4% at high densities and by only 1% at low 
densities, well within the errors. Our results will not be 
changed by choosing time dependent mass limits, at least 
for levels that correspond to any likely evolution in the 
mean mass buildup in clusters. 

4.3. Field Galaxy Infall 

Infall of late-type galaxies from the low-density region 
surrounding the cluster is expected to play a continuing 
role in the build-up of early-type cluster galaxies through 
their transformation to elliptical and SO galaxies (and 
thus leads to a build-up in the total mass in early-type 
galaxies with time). Since the early- type fraction in the 
field above our ma ss limit IO^^'^Mq is 30% lower than in 
the c luster (2±50% lNoeske et al.l[2007[ : Ivan der Wei et al.l 
|2007( ). the infalling spiral and late-type galaxies must 
transform to early-type galaxies if the infall rate adds 
significant mass to the cluster (even if the mechanisms of 
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Fig. 9.— Mass-density relation for MS 1054 and RX J0152 at 2 = 
0.83 for a magnitude-limited sample. The red circles are elliptical 
galaxies, the orange squares are SO's and the blue spirals are late- 
type galaxies. Each galaxy has My < + 1, the criterion we 
use for luminosity-selected samples. The horizontal green dashed 
line shows our mass selection limit while the vertical black dashed 
line shows the division between the two density bins in Figs. [4] and 
[7] Black circle s highlight red g alaxies that are pote ntial mergers 
in M S 1054 jTran et al.ll2005al) and in RX J0152 UPostman et al.l 
[2OOI) . The mean mass of a late-type galaxy, in our luminosity- 
selected sample, is at the green line, or IO^'^'^^'^'^Mq, one quarter 
of the mean mass, lO^^'^^^ ^M©, of the early-type galaxies in our 
luminosity-selected sample. A significant fraction of the late-type 
star-forming galaxies are below our mass limit, yet they would 
be included in the luminosity-selected, MDR for these z = 0.83 
clusters. At lower redshifts, as star formation ends and they fade 
to lower luminosities, they will drop out of the luminosity-selected 
MDR. The result would be an apparent increase in the early-type 
fraction in the luminosity-selected samples. 

the transformation are unclear at this time). Otherwise 
the fraction of early-type galaxies would decrease with 
time, which is inconsistent with what is seen in either 
the mass- or luminosity-selected samples. One caveat 
on this argument is that the low-density regions around 
rich clusters are likely to be unlike the "true" field (see 
e.g. iDiaferio et al.l[200ll ). and so studies of the cluster 
environment within ~ 10 Mpc (comoving) will be needed 
for quantitative discussions of the role of infall. 

We expect clusters of galaxies to be continually grow- 
ing in mass by the infall of galaxies and groups of galaxies 
from larger radii. Yet we find no change in the mor- 
phological mix above our mass limit. No change in the 
fraction of early-type galaxies does not preclude signif- 
icant mass growth, however. A simple argument based 
on the rough constancy of the observed ~85% early-type 
fraction allows us to estimate how the mass in early- 
type cluster members could grow. First, we assume that 
the infall rate is always high enough to keep the frac- 
tion of late-type galaxies ~15%. Second, we assume that 
galaxies take ~ 1.5 Gyr to transform from a late- type 
spiral into an early-type ga laxy, the assumption from 
Ivan Dokkum fc Franxl (|2001f ). Between z = 0.83 and 



= 0, there are 7 Gyr, slightly more than 4 times the 
transformation timescale. Every 1.5 Gyr, the mass in 
early-type galaxies in this scenario increases by 18%, 
as 15% of the total population of galaxies transforms 
into same type as the remaining 85% of the population. 
These assumptions lead to a doubling of the z ~ 0.8 
mass by 2; = 0. In other words, ~50% of the z ~ 
early-type galaxy population enters the cluster as late- 
type galaxies between z ~ 0.8 and z ~ 0. Interestingly, 
this rate of increase in the mass density of ellipticals and 
SO galaxies is roughly comparable to that typically mea- 
sured in field surv eys (see iBell et al.l l2004t iFabeil 120071 : 
iBrown et al.ll2007l for example). The time-scale of trans- 
formation and the infall rate are the largest uncertain- 
ties. A longer transformation time-scale would in turn 
lower the amount of inferred evolution among early-type 
galaxies above our mass-threshold, and this calculation 
neglects that many infalling massive galaxies would al- 
ready be elliptical and SO galaxies. 

4.4. The Role of Low -Mass Galaxies - Transformation 
at Low Masses 

Our mass-selected sample reaches down to masses 
10^°-^ Mq, and the lack of evolution of the early- type 
fraction in that sample leads us to infer that changes at 
lower masses are playing a role in the evolution seen in 
luminosity-selected sam ples, a possibility discusssed in 
iDe Propris et al.l ()2003[ ) . No comprehensive studies have 
been carried out that are complete from high to very low 
galaxy masses. 

However, there is some evidence of evolution at low 
masses (transformation from late-type to early-type) 
but the results have not been consistent. For exam- 
ple, a number of authors have arg ued for evolution in 
the faint, red clu s ter po pulation (Dc Lucia et al. '200^ 
Goto et al.l 120051: iTanaka et al. 2005; Dc Lucia et al 



2007t IStott' et al. 1 120071) . This evolution is observed as 
an increase in the ratio of faint - often called "dwarf" - 
red-sequence galaxi es to the "gia n t" ga laxies with lumi- 
nosities of L*. In iHolden et al.l ()2006f ). we speculated 
that these additional "dwarf" galaxies would become 
low-mass SO galaxies in z = clusters. Th e samples of 
IDe Lucia etall (|2007D and lStott et~al] (|2007D are selected 
to contain galaxies with masses below lO^^'^M©, so there 
is no conflict between the observed evolution in that 
sample for galaxies with masses of ~ 0.25A^* and the 
lack of evolution in our mass-selected sample containing 
> 0.25 7W*. However, the reported results are not unani- 
mous. lAndreonI (|2006D finds no change in t he population 
of faint galax ies in the c luster MS 1054. ^Tanak a et al.l 
(2005) and Tan aka et al.l (j2007f ) finds a similar result, i.e., 
no change in the luminosity function for the red-sequence 
in RX J 0152 and RDCS 125 2 9-29 27 at z = 1.24. Inter- 
estingly, ITanaka et all (|2005l . 120071) find a truncation in 
the red sequence for the lowe r mass groups of galaxies 
outside of the cluster cores. ITanaka et afl §005) and 
lAndreonI (|2006l ) both imply that one might find no evo- 
lution in the early-type galaxy fraction at even lower 
masses for the clusters in our sample, a result that can 
be verified with deeper spectroscopic samples of mass- 
selected galaxies. 

The inconsistencies in these results suggest that robust, 
mass-selected samples that are complete to substantially 
below > 10^"'®Mo are needed to establish the nature of 
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the evolution of the star-forming galaxies in the cluster 
environment. 

5. SUMMARY 

We selected a sample of clusters spanning a range in 
redshift from z = 0.023 to z = 0.83. Each cluster we 
included has approximately the same X-ray luminosity 
and velocity dispersion to mitigate the observed correla- 
tions between cluster mass estima tes and the fraction o f 
ellip tical and SO galax ies seen in Postm an et al.l ()2005l ) 
and lDesai et al.l (|2007f) . Using high-quality photometry, 
largely from iJST, we computed rest-frame optical colors 
for the spectroscopically confirmed members and derived 
mass-to-light ratios(A//i). From M/L, we determined 
stellar masses of the 674 galaxies in our sample. We find 
that these stellar masses agree well with those from dy- 
namical measurements, with a scatter ranging from 0.11 
dex at z ~ to 0.20 dex at our highest redshifts. 

All of the galaxies in our sample have morphological 
classifications. W hile the classifications come frorn a va- 
riety of sources (,Dressleijll980al: iFabrican t et al.l [2OOOI: 
iTran et al.ll2003l . l2005bl : iPostman et al.l l2d05l. aU of the 
z > 0.1 catalogs were const ructed to be con sistent with 
the classification scheme of iDressleil (jlQSOal ). We used 
these catalogs to derive the morphology-density relations 
(MDRs) in the five clusters, and to establish the evolu- 
tion in that relation for luminosity-selected samples (for 
luminosities < My -|- 1 . ) As shown in Figures [3] and [H 
we found good a gree ment for both the z = M DR from 
iDressleii (|l980l7) and IPostman fc Gelleii (1984), as well 
as for the evolution seen in previous work (.Dressier et al.l 
fl997LlSmith et al.ll2005l : [Postman et al.ll2005n . for the 452 
galaxies brighter than < My -I- 1 in our sample. 

We contrast the results for our luminosity-selected 
samples with a similar large completeness-corrected 
mass-selected samples of 441 galaxies with masses > 
VS^^-^Mq (0.25 at z = 0) in each of the clusters. 
At this mass limit, the majority (~70%) of the stellar 
mass in galaxies in the clusters is in our mass-selected 
sample. 

Very interesting differences are seen between the mass- 
and luminosity-selected MDRs. First, the MDRs derived 
for the individual clusters from the mass-selected galaxy 
samples typically have larger early-type galaxy fractions, 
by 10%-20%, and show marginal evidence of a lower 
slope with density. Second, and particularly striking, is 
that essentially no evolution is seen in the mass-selected 
early- type galaxy fraction with redshift. The change of 
the slope is substantially less than for the luminosity- 
selected sample, and is consistent with no change with 
redshift. We find mean early-type fractions, 83 ± 3% 
at 500 galaxies Mpc~^ > S > 50 galaxies Mpc~^ and 
92 ± 4% at E > 500 galaxies Mpc'^ This indicates 
that the evolution in the early-type fraction observed in 
luminosity-selected samples of galaxies arises predomi- 
nately from the inclusion of lower mass, < VS^'^'^Mq^ 
late- type galaxies at higher redshifts. 

Despite finding no evolution in the mass-selected early- 
type fraction with redshift, we recognize that the cluster 
galaxy formation process is not finished at z = 0.83. 
The overall mass in early-type galaxies is expected to 
increase, but this increase must happen so as to leave 
the early-type galaxy fraction largely unchanged. The 
early-type fraction in the field above the same mass 



limit is very significantl y lower ( 48% ± 7% at z = 0.8 
Ivan der Wei et al1l2007f ). Infall of field galaxies, with a 
larger fraction of spiral and late-type galaxies, into clus- 
ters is expected and must be accompanied by transfor- 
mation to early-type galaxies to ensure that the fraction 
of early- types remains constant. 

Mass-selected samples to low masses, along with de- 
tailed morphologies, are realizable and would be of great 
value for characterizing and constraining galaxy evolu- 
tion over the last 7-8 Gyr out to z ^ 1. Identifying 
these physical processes and the nature of the transi- 
tion from the field early-type fraction is the next step - 
for a range of cluster m asses (see e.g., IPostman et al.l 
I2005t iDesai et al.l I2007D . Of particular value will be 
the outer regions, where future cluster members will 
be transformed during their infall into the cluster by 
z ~ 0. Significant insights into the build-up of galax- 
ies in dense environments over the last 7-8 Gyr would 
result from increasing the sample of clusters with mass- 
selected galaxies that have redshifts, morphologies and 
excellent photometry. These data would be of partic- 
ular value for developing insight into the formation of 
SO galaxies which, show strong evolution in luminosity- 
selected samples from z ~ 1 to ~ 0. 

In summary, we have derived luminosities, masses and 
morphologies for a large and complete sample of clus- 
ter galaxies ranging from z = 0.023 to 0.83. The 
luminosity-selected galaxy samples show, in all clusters, 
a clear MDR. Furthermore, the early-type galaxy frac- 
tion evolves substantially with redshift, as found by oth- 
ers. When we construct mass-selected samples, the MDR 
are offset towards larger early-type fractions and are 
marginally weaker. More dramatically, the evolution in 
the early-type galaxy fraction is very weak, being con- 
sistent with no evolution over the 7 Gyr period. As 
these two samples are drawn from the same parent sam- 
ple of five massive clusters of galaxies, we conclude that 
the galaxies that lead to the evolution of the luminosity- 
selected MDR must be galaxies less than our mass-limit, 
namely those with masses < IO^^'^Mq. 
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APPENDIX 
FIRES DATA FOR MS 1054 

We used additional photom etry for MS 1054 from the Faint Infrared Extragalactic Survey (FIRES), discussed in 
iForster Schreiber et al] ()2006[) . The FIRES team used FORSl on the Very Large Telescope (VLT) to obtain the UBV 
data and ISAAC to obtain the JgHKs data. The ISAAC imaging data consists of four overlapping pointings with 77 
hr of total exposure covering a 5.'5 by 5.'3 field of view. 

The original FIRES catalog was created using a combination of the VLT data and older WFPC2 imaging in F606W 
and F814W. We elected to replace the WFPC2 data with our ACS imaging, as our data should be of higher quaUty 
and include an additional filter, zsso- To obtain accurate colors, we smoothed the ACS imaging, along with the rest of 
the VLT imaging, to the same effective seeing as the worst of the ground based imaging. We then re-binned the ACS 
data to the same scale, 0'.'0996 pixel"^, and orientation as the VLT Ks image. We used the un-smoothed Ks image as a 
detection image. All magnitudes for the colors were measured in fixed apertures on the smoothed data. The apertures 
ranged in size from I'.'O to 2'.'4 in size in increments of 0'.'2. We selec ted the aperture in whic h the signal-to-noise ratio 
for the color measurements was maximized. We used SExtractor (jBertin fc Arnout3ll996l ) for object detection and 
photometry for the FIRES imaging. 

As we have re-measured the magnitudes in all the images, we compared our resulting magnitudes with those as 
measured by the FIRES team. For the near infrared data, the scatter was less than 0.01 mag. For the B and V data, 
the scatter was around 0.02 mag, with the U showing 0.04 mag in scatter. Given that the images are the same, this 
scatter came from different choices of parameters used for object detection parameters, resulting in slightly different 
object positions and sizes. We examined the scatter for those objects labeled by SExtractor as "biased", i.e., objects 
with close neighbors in the FIRES catalog. We found that the scatter increased (at times it doubled) only for those 
objects flagged as "biased" . Thus, the difference in detection parameters is most likely what causes the scatter in 
the photometry between our measurements and those of the FIRES team. In addition to the scatter, we found a 
small offset between our magnitudes and those from FIRES. This result came from the different point spread function 
smoothing kernels used by us and by the FIRES team. This offset is less than 0.01 mag for all filters. 

Our photometric catalog was then merged with the ACS imaging catalog. We selected galaxies from the FIRES 
catalog, then matched those to identificatio ns from the ACS catal og. This catalog contains morphological types from 
([Postman et al.l [20051 ) and the colors from iBlakeslee et al.l (|2006l ). When two objects from the ACS imaging were 
identified as a single object in the FIRES imaging, we rem oved them. The final catalog contains 1221 galaxies, with 
a subset of 482 that have morphological classifications from iPostman et al.l ()2005[ ). 

We found a small 0.0 18 ± 0.012 ma g offset between our zgso magnitudes and those from'Blakesl ee et al.l (j^OOG) . The 
magnitudes in .Blakcsl ee et al.l (|2006[ ) have a 0.2 mag aperture correction to estimate total magnitudes. The FIRES 
magnitudes are based on the five Kron radius aperture magnitudes from SExtractor as measured on the Kg imaging. 
The Z850 — Ks color is then used to compute the total fiux in the zgso filter. However, no aperture correction is made, 
so this 0.018 mag offset reflects how much light is lost in the Kg imaging. 

MASS ESTIMATES BASED ON MULTI-BAND PHOTOMETRY 

For MS 1054, we estimated stellar masses by fltting spectral energy distributions (SEDs) to the photometric data 
from the FIRES data set. Specifically, we used two sets of BC03 models: simple stellar populations (SSPs) and models 
with exponentially decaying star formation rates. These were fit to the combination of the FIRES data and ACS 
imaging. For both sets of BC03 models, we assumed a Salpeter initial mass functions. When using models with 
exponentially decaying star formation rates, or r models, we assumed timescales (r) ranging from 100 Myr to 5 Gyr. 
For both the SSPs and the r models, we used three different metallicity values, 0.4, 1, and 2.5 Zq. 

Given a SED and the photometry of a galaxy in one of the clusters in our sample, we redshifted the SED to the 
observed value. We computed the fluxes in the observed filters usin g the filter curves pr ovided by ESO for FORS and 
ISAAC, and by the Space Telescope Science Institute for the ACS TSir ianni et al.l 120051 ). For a given metallicity and 
r, we allowed the range of ages to vary. For every model, we computed the normalization n that minimized the 
for that model. This eliminates the need to fit for the best model and the best normalization, as the normalization is 
the one that automatically minimizes the for a given model. The normalization is n = {J2 fi''^i / '^f) / iJ2 / '^i) 
where i denotes a filter, fi is the observed flux, rrii is the model flux and at is the estimated error. The error on n is 
Cn = (X) X] ™?)^^^- This normalization provides the mass in stars required to reproduce the observed flux in the 
color aperture. We used the total magnitudes in the Kg filter to correct this normalization to a "total" stellar mass. 
We then compared the x^ for each model using this normalization while varying the other parameters of interest. 

These model fits were also used to derive rest-frame colors. We computed the model's z = 0, or rest-frame, 
magnitudes in the desired filters, such as B and V. We then computed the model magnitudes in the observed filters 
as part of the x^ fitting. We used the offset between the best-fltting model's observed and rest-frame magnitudes to 
estimate the rest-frame magnitudes from the data for each galaxy. 

For both the SSP and r models, we included the effect of dust by adding an additional foreground screen to our 
spectral energy distribution model flts. The typical change in the resulting stellar mass was minimal. For example, the 
average mass changes by 7% between r models with dust and those without. The models generally flt the data with a 
minimal amount of dust for most galaxies, except for the bluer, star-forming objects. Removing the fllters that sample 
the rest-frame UV, the U and B in the FIRES data sample <2500Aat z — 0.83, reduces this change in the mass to 
4%. When combined with r models, dust removes any mass-metallicity relation and, as expected, moves the average 
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age to younger values. The mode of the age distribution, however, stays the same. The modal age of the population is 
the age of the "red-sequence", and those galaxies are all fit with little or no dust. The addition of dust adds a tail of 
younger ages. Thus dust models significantly changed the ages and masses only for the bluer populations, as expected 
f or galaxies with o n going star formation, and do not affect the majority of the cluster population. 

'Bel l et all JlOOf) give, in §2 of their Appendix A2 (Table 7), a re lation between rest-frame color and M/L in the 
Johnson B, V, and R passbands. We found, as in iBell et al.l (|2003f l. the best-fitting SED M/Lb correlates strongly 
with the B — V rest-frame color for the galaxies in MS 1054, see Figure [TOl This plot is made by fitting dust-free 
T models to all of the FIRES and ACS filters. We show our best-fitting rel ation betw e en M /Lb and the rest-frame 
B — V color. We find that our relations are very similar to the results from I Bell et al.l (|2003f ). shown as dashed lines 
in Figure [TUl 

We found a scatter of 33%, or 0.14 dex, in M/L around the best-fitting relation for MS 1054. Early-type galaxies 
have a slightly larger scatter of 34% or 0.15 dex, while late-type galaxies have 30% or 0.13 dex. In our lower redshift 
sample created using the SDSS DR4, we find a smaller scatter, 28% or 0.12 dex, with the same scatter for both late 
and early-type galaxies. 
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Fig. 10. — M/Lb as a function of rest-frame B — V color for the MS 1054 sample {top) and the SDSS DR4 sample (bottom). In each 
panel, blue spirals are late-type galaxies while the orange squares are early-type galaxies. The solid lines repres ent the best fit t o the 
relation between the B — V color and the mass-to-light ratio in the B band. The dashed lines are the relations from IBell eFaTI 12003) . We 
reproduce the results of lBell et al.l I I2003I) . and thus used that relation to estimate M/Lb for the whole of our cluster sample. The scatter 
around the relation is 0.14 dex in the high-redshift sample and 0.12 dex for the low-redshift sample, with no significant difference in the 
scatter for early or late-type galaxies. 



